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Non-deterministic chaos is a new dynamical paradigm where a non-deterministic system is influ- 
enced by random perturbations to produce the appearance of complexity. The non-determinism is 
envisioned to occur only at a single point in phase space, where many trajectories intersect. In the 
presence of external random perturbations (noise), whenever the phase space trajectory approaches 
, the singularity, it will jump in an unpredictable way to a different solution. This behavior, while 

0^ ' similar in appearance to deterministic chaos, has rather different implications for prediction and 

0^ ' control. 
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QO ■ INTRODUCTION 

It is often assumed that classical physical systems are deterministic , meaning that the forward time evolution is 
^ uniquely determined by the current state of the system. Correspondingly, determinism also implies that the system's 

' current state uniquely determines its past behavior. These qualities follow directly from the Existence and Unique- 
ness Theorem of differential equations which requires that the equations of motion of the system are Lipschitz 
00 continuous. However, this is nothing in classical mechanics that requires Lipschitz continuity. Indeed, in the case of a 
\ cracking whip, the physical solutions imply violation of the Lipschitz condition Q. A similar effect is seen in seismic 

■ waves as they approach the surface of the Earth If Nature is not required to be Lipschitzian, we must ask if it 
required to be deterministic. Some types of non-determinism in classical systems have been explored previously [^^. 

■ In this work, we examine a somewhat different flavor of non-determinism, one which implies the possible existence of 
a new dynamical behavior: non- deterministic chaos. 
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^ . THE NON-DETERMINISTIC HARMONIC OSCILLATOR (NDHO) 

O 
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We begin with the simple harmonic oscillator (SHO), described by the equations 



y 



C3 ■ dt 

Jt" ^ 

Solutions of eqs. |l| are circles in the (x, y) phase plane. The SHO is a deterministic system, so that every point (x, y) 
belongs to a unique solution described by a circle with a particular radius r — \/ x^ + y^- 

Suppose, now, that we apply the following non-linear coordinate transformation to the SHO phase space: 



a; — > X — r = a; 



V^^ + y^- (2) 



This translates all points on a circle of radius r in the positive x-direction by an amount equal to r. A family of 
circles concentric about the origin in the original space will now share a common tangent point at the origin of the 
transformed space (see Figure |l|) . 

The equation of a transformed circle in the new space is given by 

(x-r)2+y2^r2, (3) 
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or, solving for r, 



r^^{x^+y\ (4) 
2x 

Using eqn. ^, we can easily apply the transformation to the SHO. The transformed SHO equations of motion in the 
new coordinate system are given by 

d 

dty-2x 2"^- 

From the above discussion, the solutions of eqs. |^ will be the family of transformed circles, all sharing a common 
tangent point at the origin. Such intersection of many phase space trajectories is not so unusual. An attracting fixed 
point, for example, is approached asymptotically for all initial conditions in its basin of attraction (i.e., solutions are 
unique for finite times). What is unusual about eqs. ^ is that the common point is intersected in finite time, and 
further is not a fixed point. This is easily seen by taking the limit of eqs. ^ along a solution of radius r: 

d „ 
lim —X = 0, 

x,y^a dt 

lim —y = r. (6) 

x,y^o dt 

Thus, the origin is a singularity of eqs. |^, where neither past nor future time evolution is uniquely determined. Hence- 
forth, we shall refer to eqs. |^ as the non- deterministic harmonic oscillator (NDHO). 

The NDHO provides the paradigm for the type of system we are examining: solutions of the equations of mo- 
tion are a family of closed loops ("transients") all sharing a common tangent point. From eqs. ^ the dynamics of 
the NDHO are not defined by the equations of motion alone (this is not a necessary condition for a system to be 
non-deterministic Q). However, let us imagine we built an NDHO in a laboratory. How would it behave? We note 
that all physical systems are subject to external perturbations, or "noise" . While the physical state of our NDHO 
is far (in phase space) from the point (0,0), external noise will have little effect, provided the average amplitude of 
the fluctuations is small compared to r for that trajectory. However, as the trajectory approaches the origin, noise 
plays a larger role. Solutions for all r converge together, ultimately intersecting at (0,0). Thus, noise will cause the 
trajectory to jump between solutions of widely differing r in a random way. 



What is the effect of this on our laboratory NDHO? Suppose we begin the system on a solution of radius ri. As we 
watch the system evolve forward in time, we will find that after it passes near the origin the trajectory has changed 
to a completely different solution of radius r2. Repeating the experiment with the same initial conditions, we find 
that the trajectory jumps to a completely different solution of radius r^, where r^ ^ r-i- Were we to repeat this a 
large number of times, for different values of ri, we would find that the solution after the singularity is completely 
unrelated to the solution before. If the NDHO were allowed to run for several oscillations, a time series measurement 
of one variable would appear as a piecewise continuous sequence of oscillations with different amplitudes. Further, 
the sequence of amplitudes would be random and unpredictable. We term this behavior non- deterministic chaos, 
non-deterministic because its origin lies in the non-determinism at a non-Lipschitz singularity, and chaos because of 
the long term unpredictability of the dynamics. 



A PHYSICALLY MOTIVATED EXAMPLE 



The NDHO, while illustrative of the type of non-determinism we are examining, is also a somewhat contrived ex- 
ample. We now describe a non-deterministic system based on physical considerations. This system is a model of the 
behavior of neutron star magnetic fields. We describe it briefly; for a more detailed discussion, the reader is referred 
to §. 

The model envisions two oppositely charged spherical shells which are allowed to rotate differentially. The mag- 
netization of one shell (Mi) will interact with the magnetic field of the second shell (H2), as well as experience 



2 



non-electromagnetic ("mechanical") interactions with the surrounding medium. The magnetic interactions include a 
term to induce precession of Mi about the instantaneous direction of H2, and the Landau-Lifshitz magnetic damp- 
ing, which tends to align Mi with the direction of H2. The mechanical interaction is taken as a simple damping, 
proportional to the difference in angular velocities of the two spheres. Parameterizing the interactions, we obtain the 
following equations: 

^Mi = 7(Mi X H2) - ~X{^\-^Mi - H2) - 7f • K - UJ2). (7) 
at M( 

Following the scaling procedure described in and conserving angular momentum, we arrive at the following 
equation: 

d ^ f m. ■ z \ - ^ 

— m — — m X z — A r-m — z — e(m — z), (8) 

ar \ / 

where m is the scaled magnetization, r is the scaled time, A is the scaled Landau damping parameter, and e is the 
scaled viscous damping parameter tensor. 

Examination of eqn. ^ reveals axial symmetry about the z-axis. This prompts us to make the following transfor- 
mation: 



which implies 



X - 

z = mz (9) 
= arctan — - 



rrix X cos (p 

my — s- X sin </) (10) 



Substituting the above transformations into eqn. H, we obtain 

Xxz 



x"^ + z^ 

■ - 

x"^ + z"^ 

and 



iz - (A - -e) (11) 



0=1 (12) 

where an overdot again represents differentiation with respect to the scaled time r. The cf) equation is trivial, simply 
representing a constant precession about the z-axis. Any interesting dynamical behavior must occur in eqs. |ll|. 

Numerical integration of eqs. |ll| for e < A yields the phase space plot in Figure |[ and the time series in Figure ||. 
Note the apparent intersection of trajectories at the origin, indicating possible non-determinism at that point. Indeed, 
it can be shown that for the RHS of eqs. the point (0, 0) is a non-Lipschitz singularity. Given this, the non- 
deterministic nature of eqs. |ll| at (0, 0) is proved by the following argument (for more details the reader is referred 
to |,|): 



1. The only fixed point of eqs. 11 for e < A is at (0, 1) 



2. The sole fixed point is a saddle, and thus is not an attractor for a set of initial conditions of non-zero measure [ pl)| . 
This also implies that no periodic orbits exist pO| . 

3. By the Poincare-Bendixson Theorem, if the solutions of a two-degree of freedom dynamical system contain only 
Lipschitz points, then the only possible bounded asymptotic solutions are stationary or periodic. As neither 
possibility exists for eqs. |ll|, (0, 0) must be contained by all bounded solutions, and therefore this point is 
non-deterministic. 
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NON-DETERMINISM AND PREDICTABILITY 



Non-deterministic chaos has the property of being predictable for short times (between intersections of the singu- 
larity), yet completely unpredictable over long time periods. Long-term unpredicability is also one of the hallmarks 
of deterministic chaos, but it is here that the similarity ends. Aside from being described by deterministic equations, 
deterministic chaos is often characterized by exponential divergence of initially close solutions, and associated with 
an complex fractal structure, the strange attractor. Non-deterministic chaos derives its unpredictability from a more 
violent, but localized instability. Further, there exists no attractor, strange or otherwise, at least in the usual sense 
of the word. 

Figures H] illustrate the loss of information at the singularity due to the presence of external noise. A set of initial 
conditions in evolved via eqs. |ll|. Initially, we see the phase space stretching, similar to what occurs in deterministic 
chaos. However, once the singularity is encountered, points are scattered, and soon are randomly spread through a 
region of phase space (in this case, the region is enclosed by the homoclines of the saddle at (0,1) 0). This behavior 
is in stark contrast to what one expects from deterministic chaos, where an initial volume is stretched and folded [0, 
spreading across the attractor in a smooth fashion. Dynamical measures such as the Lyapunov exponent are meaning- 
less. The "attractor" exists only in a statistical sense, representing the probability density that a particular point in 
phase space will be visited. Formally, one could find this probability density by transforming eqs. |l^ into a Brownian 
motion and numerically solving the forward Kolmogorov equation |L2| . This approach is complicated by the existence 



of the singularity. Instead, we simply integrated eqs. [lljfor 200 million time steps, and constructed the distribution 
from this. The results are shown in Figure ^. 

Figure ^ gives a long-term, global statistical picture. However, the nature of non-deterministic chaos allows us to 
easily extract more useful statistical information. In particular, we shall utilize the fact that away from the singular 
point, the dynamics is quite well-behaved. Let us return to the NDHO, as its solutions are known analytically. The 
solutions of the NDHO may be parameterized by their radius. Solutions away from the singularity have essentially 
constant radii; the big jumps occur only near the singularity. Can we predict the probability that a circle of given 
r is chosen when the orbit leaves some neighborhood about the origin? Let us define this neighborhood as a disk of 
radius S, and note that an orbit leaving this neighborhood does so with angle 6, which we take as measured from the 
y-axis. Now, assuming the external fiuctuations to be isotropic, the probability density of picking a particular 9 is 
constant, i.e., 

p{e)de oc de. (13) 

Next, we note that everywhere except at the origin the Existence and Uniqueness Theorem applies to solutions of 
eqs. (P), thus each circle of radius r is associated with a unique 9, and we may write 6* as a function of r. Substituting 
into p(d) , we find 

d9(r) 

p{r)dr oc ^ dr. (14) 
or 

The probability of getting a circle between r and r + Ar is simply 

P{r,Ar)oc -^dr ^ 9{r + Ar) - 9{r). (15) 



dr 



For the case at hand, we find 
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Plr.Ar) oc arccos -7 — r — arccos— . (16) 

^ ' ^ 2(r + Ar) 2r ^ ' 

This approach (first described in ||^) is somewhat simplified. A rigorous derivation would account for the statistical 
properties of the noise, and derive Ar) via stochastic calculus. The above does show, however, that the sim- 
ple structure of the solutions of a non-deterministic system lends itself to the construction of statistical arguments. 
Further, with a judicious choice of (5, based on knowledge of the average amplitude of the fluctuations, the above 
procedure should yield a good approximation of the true distribution P(r, Ar). 
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CONTROLLING NON-DETERMINISTIC CHAOS 



The control of deterministic chaotic systems using smah perturbations has been a subject of recent vigorous re- 
search The most popular method of controlling deterministic chaos involves the stabilization of (otherwise) 
unstable periodic orbits which are embedded in the chaotic motion. As there exist an infinity of orbits, a rich variety 
of behaviors may be extracted from the controlled deterministic chaotic system, allowing for flexibility and easy op- 
timization of a system's behavior. 

For a non-deterministic chaotic system, we have a similar situation. With a continuum of different solutions inter- 
secting at a single point, we can easily effect control via an appropriate perturbation. Similar to the previous section, 
we simply examine how solutions leave a (5— neighborhood about the singularity. Again, away from the singular point 
the solution is well-defined. Suppose that each different solution may be parameterized by some quantity 7 (in the 
case of the NDHO, this is the radius). A given solution, parameterized by 70, will intersect the ^-neighborhood at a 
unique point {xo,yo)- From this, we may construct the angle ^(70) = arctanxo/j/o- 

The angle ^(7) we term the control angle, and the reason should be obvious. To keep the system on a solution with 
parameter 70, we need only to wait until the trajectory approaches the origin, and then perturb it so that it leaves at 
angle 0(70). This perturbation will be quite small, of the order of S. with the size of S being determined largely by 
the noise amplitude. We see that in a non-deterministic system, there is a continuum of possibilities available through 
small control perturbations. If a change in system behavior were required, it is easily and quickly effected by simply 
changing ^(7). In fact, one could vary 0(7) as a function of time to induce arbitrarily complex behavior. 

As an example, we have applied this control algorithm to eqs. |ll|. To simulate the effect of noise, a small (10~^ 
of the integration stepsize) normally distributed random number was added at each integration step. The controlled 
signals for various values of 9 are shown in Figure 0. Figure M shows the effect of noise for different values of 6. 



DISCUSSION 



The type of "non-determinism" described above should not be construed as implying stochasticity. Indeed, the 
behavior of both the NDHO and neutron star model are uniquely determined away from the singular point. It is at 
this point, and this point only, that the non-deterministic nature of the equations arises. In the presence of random 
fiuctuations, which are ubiquitous (though perhaps small) in physical systems, the non-determinism, albeit it at a 
single point, becomes important. The resulting dynamics, which we have termed non- deterministic chaos, consist of 
a random sequence of "transient" oscillations. 

This work does not represent the first suggestion that non-determinism exists in classical systems. The terminal 
dynamics described by M. Zak |^ utilizes a similar mechanism where multiple trajectories intersect at a common equi- 
librium point in finite time. Chen [|| has independently suggested the same behavior under the heading of noise induced 
instability (though non-determinism as such is never explicitly mentioned). For clarity, we shall henceforth refer to 
terminal dynamics as "Type I non-determinism" , and the dynamics described herein as "Type II non-determinism" . 
The primary difference between the two is that for Type I non-determinism, the singularity occurs at an equilibrium 
point of the equations of motion, while for Type II, the singularity is shared among a group of dynamic trajectories. 
The physical implications of this difference are yet to be explored. Other aspects of non-deterministic systems, espe- 
cially in the presence of noise, have been explored by A. Hiibler l]!^ . 

There is no principle in nature that precludes the existence of the types of systems we have described. Indeed, 
an analog VLSI circuit has been built which displays Type I non-determinism p5| ] . The neutron star model displays 
Type II non-determinism, and is based on perfectly reasonable physical assumptions. Work is in progress to build 
an electronic analog of this model. It has also been suggested that non-determinism may play an important role in 
biological systems |Lq |. However, the question arises as to the ubiquity of such systems. Is non-determinism a generic 
property of Nature, or have we simply stumbled upon a few pathological examples? We cannot at this time answer 
that question. 

However, we feel further investigation is warranted. In particular, application of several standard measures (power 
spectrum, Lyapunov exponent, etc.) to a time series generated by eqs. |ll| would lead one to believe that one examining 
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an instance of deterministic chaos |j. As we have seen, though, issues of prediction and control would be addressed 
much differently for a non-deterministic chaotic system. Indeed, Crutchfield has shown that in the context of model 
building, assuming determinism when the underlying process is non-deterministic leads to undue complexity in the 
model |17[] . It would seem reasonable to search for non-deterministic chaos in apparently complex systems, especially 
in cases where traditional analysis tools (which, again, assume determinism) have failed. In a forthcoming paper, we 
will address the problem of detecting non-determinism in observed data. 
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FIG. 1. Some examples of circular orbits of different radii, all sharing a common point at the origin. 



FIG. 2. Phase plot of solutions in the neutron star model for A — 1.0, e — 0.6, e — 0.7. 



FIG. 3. Time series z{t) vs. t for the neutron star model with parameter values A — 1.0, e — 0.6, e = 0.7. 

FIG. 4. Loss of information due to the singularity, a) 10,000 initial points are arranged in a 100x100 square, b) Initial 
evolution, c) When the singularity is encountered, initally close points are scattered randomly, d) All information about the 
initial conditions has been lost. The only information carried by the system is in the density of trajectories. 

FIG. 5. Probability of finding the system in a particular region of phase space. The distribution was found by integrating 
the equations for 200 million steps and totalling the amount of time spent in each region. 



FIG. 6. Examples of output from the neutron star model when the control algorithm is applied. Signals are shown for ^(7) = 
a) 0.005, b) 0.03, c) 0.08, and d) 0.2. 



FIG. 7. The control algorithm begins to break down if 5 is chosen to be comparable to the noise level. Signals are shown for 
(5 = a) W^a, b) lO^CT, c) lO^cr, and d) lOcr, where a is the RMS of the noise. 



6 



